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Potential formulation of the dispersion relation for a uniform, magnetized 
plasma with stationary ions in terms of a vector phasor 

Robert W. Johnsor0 
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The derivation of the hehcon dispersion relation for a uniform plasma with stationary ions subject to a 
constant background magnetic field is reexamined in terms of the potential formulation of electrodynamics. 
Under the same conditions considered by the standard derivation, the nonlinear self-coupling between the 
perturbed electron flow and the potential it generates is addressed. The plane wave solution for general 
propagation vector is determined for all frequencies and expressed in terms of a vector phasor. The behavior 
of the solution as described in vacuum units depends upon the ratio of conductivity to the magnitude of the 
background field. Only at low conductivity and below the cyclotron frequency can significant propagation 
occur as determined by the ratio of skin depth to wavelength. 
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I. INTRODUCTION 

In this article we reexamine the derivation of the heli- 
con dispersion relation in terms of the potential formula- 
tion of electrodynamics. Under the same approximations 
as used in the standard derivation, namely a uniform 
plasma with stationary ions subject to a constant back- 
ground magnetic field with vanishing thermal stress and 
space charge density, the linearized equation of motion 
for the electron flow describes only the electron cyclotron 
resonance when the ion contribution to the material re- 
sponse is included through the friction term. Only by ad- 
dressing the nonlinear self-coupling of the perturbed elec- 
tron flow to its own potential through the Lorentz term 
can a more interesting solution be found. The plane wave 
solution for a general propagation vector is determined, 
whose frequency depends upon its inclination from the 
plane orthogonal to the direction of the background field. 
These solutions are represented in terms of a vector of 
phase factors in addition to the oscillatory phase. Sev- 
eral cases of interest are evaluated explicitly. The main 
result is that propagation occurs only for frequencies in 
a range below the cyclotron frequency when the ratio of 
conductivity to the magnitude of the background mag- 
netic field is sufficiently low that the skin depth exceeds 
the wavelength in the plasma. 

The history of this derivationfi goes back several 
decades.—^— Chen surveys the experimental literature in 
his contribution-? to High Density Plasma Sources, and 
with Boswell reviews the development of the theory 
throughout the twentieth centuryj^i^ Recent examples 
have appeared in Physics of Plasmas of the application 
of this theory to a toroidal vessel^ and to an annularly 
bounded discharge chamber. l^-^l What these derivations 
have in common is their use of the electron equation of 



motion to determine the electric field rather than the 
electron flow as a consequence of their neglect of Gauss's 
law. In the potential formulation of electrodynamics, the 
electric and magnetic fields are recognized as auxiliary 
expressions describing the spatial and temporal varia- 
tions of the four-potential, which is determined by the in- 
homogeneous (source-bearing) Maxwell equations in con- 
junction with the gauge condition expressing the conti- 
nuity of the potential. 

This paper is organized as follows. First we will review 
the standard derivation of the "fast" and "slow" helicon 
modes. We will then reconsider the derivation in the po- 
tential formulation, including both the ion contribution 
to the material response and the nonlinear self-coupling 
of the electron flow to the potential it generates. An 
example of the use of the potential formulation for the 
electrostatic case is given by Jankauskas and KvedarasJ^ 
The solution of the material equation of motion is ex- 
pressed in terms of a vector phasor describing the direc- 
tion of the electron flow in addition to the scalar phasor 
describing its magnitude. The solution of the potential 
equation of motion is expressed in terms of a complex 
propagation vector describing the wavelength and decay 
of oscillations at a given frequency. These vectors are 
evaluated explicitly for a range of frequencies surround- 
ing the electron cyclotron resonance, and the behavior 
of the solution is found to be determined by the ratio of 
the material's conductivity to the magnitude of the back- 
ground magnetic field. We will close by discussing how 
the theory must be extended before it can be applied to 
the description of an actual experimental configuration. 



II. LINEARIZED DERIVATION 



First let us look at the derivation of the "fast" helicon 
^'robjohnson@alphawaveresearch.corii] mode, followed by the derivation of the coupled "fast" 
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2 



is based on the field equations 

V X E„ + dtB^ = , V • B„ = , 
V X = fio3^ , 

and the Unearized material equation of motion 
Ec^ = Jt^ X Bo/eno , 



(1) 
(2) 



(3) 



subject to the constraint V • J^; = 0, where the sub- 
script uj identifies the oscillating quantities, the constant 
background field Bq = BqZ defines a direction z, the 
uniform plasma density no equals the number of elec- 
trons required for neutrality, and the plasma current 
J„ = —enoV^ for electron fluid velocity V^^. The electric 
field is eliminated by the substitution 

- dtB^ = V X (J„ X Bo)/eno (4a) 
= (Bo • V)J^/eno , (4b) 

where Vno = 0, and the current is eliminated by 

-atB„ = (Bo-V)(VxB„)/^oeno, (5) 

which for a traveling wave with phase exp* (k • r — ujt) oc 
exp'^{kzZ — ujt), using the notation exp*((S) = e"^, yields 
the relation 



0= (a- Vx)B^ 



(6) 



where the total wave number a — {uj/kz){p.Qeno/BQ) 
equals the magnitude of the propagation vector a = |k|, 
leading to the Helmholtz equation (V^ + Q;^)Bt^ — de- 
scribing the "fast" helicon mode. 

For finite conductivity represented by a collision rate 
h' and including the effect of inertia represented by the 
electron mass rrie, the material equation of motion be- 
comes 

= (j^ - ia;)me/e^no + x Bo/eno , (7) 

whose curl leads to the relation 

0^{[{uj + iiy)/ik,Lj,)] (Vx)2 - (Vx) + a} B„ , (8) 

in terms of the electron cyclotron frequency ujc = 
eBo/rrie, which can be factored^ as 

=(/?+- Vx)(/?_ - Vx)B^ . (9) 

The two roots for nontrivial B^j, in terms of 7 = 
kzU!c/'2{u! + iv), give total wave numbers of 



/3± 



1 ± (1 - 2a/7) 



1/2 
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(10) 



which are identified as the "slow" and "fast" modes, also 
knownSii^ as the Trivelpiece and Gould mode for (3+ and 
the helicon mode for /3_ respectively. 

Let us now look at what Gauss's law has to say about 
these models, whose inclusion turns Eqns. ([T|) and ^ 



into the system for pre-Maxwell electrodynamics. For 
either Eqn. ([3]) or Eqn. ([7]) one can write 

V • E„ = V • (J^ X Bo)/eno (11a) 
= Bo • (V X J^)/eno , (lib) 

as the source generating the background field is external 
to the region of consideration, whereupon eliminating J^^ 
gives 

V • E„ = Bo • [(Vx)2B„]/Aioeno (12a) 
= -V2(z-B„)(So/Aioeno) . (12b) 

For consistency with the approximation of neutrality 
V • Etj = 0, one requires either fc^ — 0, where k = |k|, 
or z • Bt^ ~ 0, which then implies J^j = Ju>i yielding 
kz — 0. For Eqn. ([3]), either condition results in the ex- 
pression ijj — equivalent to the statemenl^ "that the 
magnetic field is constant [] inside a perfect conductor." 
For Eqn. ([7]), Faraday's law can be written 

- 9tB^ = (Vx)^B„(j/ - iuj)me/ fioe^no 

+dziV xB^){Bo/^ioeno) , (13) 

which can be simplified to 

iujB^ = {Bo/^Qeno)[k'^{iy - iuj)/uJc - ^zkxJB^^ . (14) 



(15a) 
(15b) 



Under the condition kz = one has the relation 
/c^ = ifiQUJUJceno/ Bo{iy — ioj) 

in terms of the AC conductivity a = e'^no/me{i' — ii^), 
which one recognizes as the usual dispersion relation for 
a conductor"'^ - up to its neglect of the term /ioco^^ arising 
from the displacement current. This result is consistent 
with the observation that the parallel current along z is 
impervious to the effect of the linearized Lorentz term 
appearing in the material equation of motion. 

The model as presented in the literature is quite spe- 
cific about its neglect of the displacement current in the 
plasma region. As Yano and Walker state, "the dis- 
placement current in [the Maxwell- Ampere equation] is 
neglected for calculation of the plasma field, as it is al- 
ways negligible in experiments." Let us consider the ef- 
fect of its inclusion. From Eqn. ([7]), the plasma current 
can be written as 



•iui = eno [{u — iuj)nie/e — Box] ■ 
-1 



= (T 



1 e 

1 
1 



•E, 



(16) 
(17) 



where ^ = uJcjiy — iw), which defines the gyrotropic con- 
ductivity tensor. 



Z • E,, 



1 
e 1 
1 + 



E„ 



(18) 
(19) 
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(20) 
(21) 

(22) 
(23) 



The Maxwell- Ampere equation can thus be written as 

V X B„ = ^o'^odt'Euj + /-ioJtd 
= Aio(eo9t + Z) • E„ , 

whose curl yields the relation 

- V^B^ = ^lo{^o^t + Z) ■ (V X E^) 
= -Hodtieodt + Z) • B„ , 

whereupon substitution for harmonic oscillations gives 

0={~k^ +uj'^/cl + ifioujY.)-B^ (24) 
= T • B^ , (25) 

for /uqCoc^ = 1. For nontrivial B^, one requires the ma- 
trix T to be singular (non-invertible), thus its determi- 
nant must vanish dot T = 0, yielding the dispersion rela- 
tion 







(26) 



which has three positive solutions indexed by 77 G 
[— 1, 0, 1] and expressed as 



= w^/cq -t- i^iQUJa/ (1 + i^r]) . 



(27) 



The solution fco is commonly identifiedi^ as the ordinary 
mode, and the solutions k± as the extraordinary modes. 

The approximation of neutrality, according to 
Eqn. (jllap . now requires 



cx (-V2 + df/cl) (z • B, 



(28) 



where the displacement term has contributed to the first 
factor. The condition B^ _L z requires (eo9t-|-5I) -E^ || z. 



thus E^ 



z and k ± z as before. The condition kz 



w^/cq can be satisfied by ct — > 0, i.e. the vacuum, or for 
the extraordinary modes k± at the frequency uj = —iv, 
which describes not an oscillation but a solution that 



decays exponentially with time, e 



Note that 



the satisfaction of Gauss's law for a neutral medium is 
what prevents these models from supporting propagation 
with a component parallel to the background field. 



III. NONLINEAR DERIVATION 

Our first objection to the preceding derivations of 
the helicon dispersion relation is mathematical, as the 
linearized form of the Lorentz term neglects the self- 
interaction between the electron fluid and the potential it 
generates. In order to describe properly the phenomenon 
of resonance within the electron fluid, one needs to ac- 
count for the influence one fluid element has on another. 
That effect appears formally as a nonlinear contribution 
to the Lorentz force which can dominate the dynamics 



when the background field is weak compared to that gen- 
erated by source currents within the plasma. 

Our next objection is not so much one of mathematics 
but rather one of physics, in that the equation chosen to 
determine the electric field, either Eqn. ([3]) or Eqn. ([7]), 
does not tell the whole story of the material response. 
The specification of charge neutrality requires the pres- 
ence of a positively charged ion background in addition 
to the negatively charged electron fluid. Just because its 
fluid velocity vanishes does not mean that its equation 
of motion is worthless, as the ion background must still 
satisfy the equation for its contribution to momentum 
conservation. 

Our last objection to the prevailing derivations is 
that the abbreviated set of Maxwell equations found in 
Eqns. ([TE]) does not satisfy the formalism of the poten- 
tial formulation of electrodynamic field theory i^^i^^ Con- 
spicuous by its absence is Gauss's law, which describes 
the fundamental relation between a source and the field it 
generates and is equivalent to the Maxwell- Ampere equa- 
tion under a Lorentz transformation. Also absent is the 
term for displacement current vital to the description of 
electromagnetic oscillations as well as the covariant ex- 
pression of the continuity of the source terms j^^^ Let us 
now reconsider the derivation of the dispersion relation in 
the potential formulation, based upon the complete clas- 
sical Maxwell system, for a model which includes both 
the contribution from the ion equation of motion as well 
as the nonlinear self-interaction of the electron fluid. 



A. Model equations 

For each species s £ {i,e}, continuity is expressed 
through the total derivative of density defined by 



hs = dtUs + V • {usYs) , 



(29) 



where Vg is the species flow velocity and fis is the particle 
source rate, and the convective derivative of velocity is 



V, = dtY, + (V, • V)V, 



(30) 



Considering a singly ionized species such that Ci — e 
and = — e under the assumption of vanishing ion 
flow Vi = 0, both the net current and net momen- 
tum are proportional to the electron flux, neVg = 
— J/e = K/me, thus the fluid acceleration can be written 
K = TOe(neVe -|- TieVe). The Statements dtVi — and 
dfjii — imply that the ions provide a fixed background 
for the electron dynamics. The equations of motio n^^i^^ 
in the absence of particle sources fis = are thus 

V • (n,T,) = F,e + en,E , (31) 
meUeVe + V • (n^Te) = Fe, - ene(E + X B) (32) 

where n^Ts is the thermal stress tensor for species s. 
Just because the ion flow vanishes does not mean that 
its equation is worthless, as the electron flow makes an 
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appearance through the friction term, F^e = —F^i = 
nienei^ei^e for interspecies coUision rate j/gi. Summing 
Eqns. dSI]) and §^ yields 

TOerieVe + V • {uiT i + UeT e) = e{ni - ne)E - erieVe X B , 

(33) 

which is simply the statement of net momentum conser- 
vation in the plasma under the given conditions, 

k + V • (nT) = jE + J X B , (34) 

where j = X^s ^s^s is the net charge density and nT ~ 
X^ HgTs is the net thermal stress. Under the approx- 
imation of uniform thermal stress V • (risTs) = and 
enforcing neutrality rig = rii such that V • (neVg) = 0, 
one is left with the system of equations 

- meVei'Ve = eE , (35) 
mAdt + Ve • V)Ve = eB x Ve , (36) 

which displays the electron fluid's coupling to inertia 
through TOe on the LHS and to electromagnetism through 
e on the RHS. Only by satisfying both equations of mo- 
tion does one achieve a consistent theory for the material 
response. 

Turning now to the electromagnetic sector, the homo- 
geneous Maxwell equations are satisfied identically in the 
potential formulation B = V x A and E = — V$ — 9(A, 
thus they cannot determine any physically relevant de- 
grees of freedom. The physical relation between the 
potential and its source is given by the inhomogeneous 
Maxwell equations 

V X B - AioeoStE - mqJ , V • E = j/eo , (37) 

where the inclusion of the displacement current is essen- 
tial to the description of electromagnetic propagation as 
well as the covariant expression of charge conservation. 
Selecting the Lorenz gauge dta + \7 ■ A = 0, where the 
scalar potential a = ^/cq, the continuity of the poten- 
tial mirrors that of the source dtj -I- V ■ J = 0, and the 
field equations are covariant —O'^A'^ — iiqJ'^, where the 
d'Alembertian operator is defined as = — Cq'^O^. 
For a neutral medium j = the scalar potential vanishes 
a = (technically is constant), thus one is left with the 
gauge condition V • A = and the field equation 

- V^A + c^^dfA = ^oJ , (38) 

expressing the relation between the vector potential A 
and its source current J. 

B. Linearized solution 

Considering a region with uniform plasma density 
Vne = such that V • Vg = and subject to a con- 
stant background magnetic field Bq along z, the lin- 
earized equation of motion reads 

dtVe = (e/me)Bo x Ve . (39) 



Nothing defined in the model breaks translational in- 
variance, so let us work in Cartesian coordinates with 
position vector r — (x,y,z), and let us write the 
phasor for the electron flow Ve = ReVe as Vg — 
eyip'^{—ujt){Vx,Vy,Vz)- The phasor equation of motion 
in terms of the cyclotron frequency ujc is then 

{V,,Vy,V,) ^ iiu,/u){~Vy,V,,0) , (40) 

thus the flow along z must vanish Vz = 0, and the re- 
lations Vx/Vy = —i(wc/uj) and Vy/Vx — iiuJc/to) must 
hold independently. For non- vanishing flow, one can then 
write 

= 1- {yx/Vy){yyiVx) (41a) 

= 1 - {L0,/L0f , (41b) 

with the positive root u — Uc which describes the mate- 
rial response at the electron cyclotron resonance. Having 
assumed away most of the relevant physics, there is noth- 
ing else left for that equation to describe. The phasor 
solution can thus be written in terms of a scalar phasor 
exp*((5) for the magnitude and a vector phasor f for the 
direction, 

Ve = Fexp*(^)f (42a) 
= Fexp*(-Wei)(l,i,0) , (42b) 

where 1^ is a real constant bearing units of velocity and 
f = -I- jf_ is the complex sum of real unit vectors sat- 
isfying the requirements of normalization |f_|_| = |f_| = 1 
and orthogonality f+ • f- =0. Note that this solution 
does not describe a propagating wave as the entire re- 
gion is oscillating at the same phase in time, having no 
dependence on position r. The vector phasor f describes 
the direction of the fiow at two times separated by one 
quarter of a cycle with a relative phase exp*(— 7r/2) = ~i. 

C. Nonlinear solution 

Let us now consider the coupled nonlinear system of 
equations 

a^Ve + (Ve • V)Ve = (e/me)(Bo -f V X A) X Vg (43) 
-V^A + c^^d'iA = -AioeneVe , (44) 

where the phasor for the generated potential A = 
{A/V) exp*((5/i)Ve has a phase of 5a relative to Ve, and 
seek solutions for Ve and A under the constraints 

V-Ve = 0, V-A = 0, (45) 
dtA - (mei/ez/e)Ve = , (46) 

where the final constraint is Eqn. psp in terms of the 
potential. The oscillatory phase 5 — k • r — cji now allows 
for variation in space through the complex propagation 
vector k = k-)_ + «k_ whose real part k-|_ determines the 
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FIG. 1. (Color online.) The propagation vector k has an 
inclination from the plane orthogonal to z of ())± determined 
by the ratio of the oscillating frequency uj to the electron 
cyclotron frequency ujc- Shown are the vectors k for cf>+ as □ 
and (f>- as when uj/ujc = 3/4. 



oscillation's wavelength and whose imaginary part k_ 
describes its attenuation. The propagation vector can 
also be written in terms of its direction, 

k = (cos Ok cos (j)k , sin Ok cos , sin cjife ) , (47) 

and magnitude k = k exp*(5fe) as k = fck, where the angle 
Ok e [— TT, tt) is its azimuth from x, and the angle (pk £ 
[— 7r/2, 7r/2] is its inclination from the plane whose normal 

is z. In these coordinates V — ?> ik, and the divergence 
constraints in Eqns. (l45t imply k • f = so that the 
convective acceleration vanishes (Ve • V)Ve = 0. 

The divergence constraints can be satisfied construc- 
tively byjipecifying = -(k^Jx + kyfy)/k^ for both 
Ve and A. The remaining constraint requires wA = 
i{mei^ei/e)'Ve, thus determining — vr/2 and A/V — 
meVei/ei-^- In this formulation, the ion contribution 
Eqn. (pS)) plays the role of Ohm's law, and the electron 
contribution Eqn. (|43p . rewritten with the nonlinear term 
isolated as 



[imje)dt - Box] V, - (V x A) x V, , 



(48) 



describes the oscillatory equilibrium in terms of momen- 
tum balance. One factor of e** appears on either side 
of that equation, but a second one appears on the RHS, 
which can be rewritten as (V x A) x Ve = 7exp'(25)k, 
where 7 = AI/[ fy]^[fx A l"^ in terms of 



r = 



1 + cos^ 9k cot^ <j)k cos 9k sin 9k cot^ (f>k 
cos 9k sin 9k cot^ 1 4- sin^ 9k cot^ (pk 



(49) 



The common factor of 7 lets one use the z component of 
Eqn. (gg), 

(me/e)5tV; =7exp^(2<5)fc, , (50) 
to simplify the x and y components into a system of 



equations for fx and fy, 

cos Ok sin^fc cot^ (f>k + i{ijJc/<^) 



fx/fy 



fv/ fx 



1 -I- cos^ 9k cot^ 4>k 
COS 9k sin 9k cot^ 4>k — i{<^c/<^) 



1 + sin^ 9k cot^ (j)k 
which one can solve for (pki^^) by writing 

0=1 -(/.//,)(/,//.) 

cx {u!c sin (pk — ^) (wc sin 0^ -t- w) , 



(51) 
(52) 



(53a) 
(53b) 



with the solutions (j)± — ± arcsin(a;/wc) as shown in 
Fig.[Tl which says that propagation at a frequency u < lUc 
is mostly along ±z while at low frequency < a; is mostly 
perpendicular to z. The ratio 



fz/ fx 



[cos 9k + i juJc/oj) sin 9k] cot (pk 
1 + sin^ 9k cot^ (pk 



(54) 



in terms of (p± such that cot(j)± — -^(yj^jiij^ — XYl"^ lets 
one express the vector phasor as 



\ + [u?Ju? -\)sxn^ 9k 
(1 — uol/uo'^) cos9k sin9k + i{ujc/(jj) 
^(uj\l(j? [cos 9k + i {ijJc/^) sin Ok] 



, (55) 



where the constant of proportionality normalizes both 
the real and imaginary components of f , which is a pri- 
mary result of this investigation. One can verify that 
k • f = so that continuity is preserved. 

Let us examine two specific cases of interest. Suppose 
first that 5i = kz — cut. The explicit solution of Eqn. p8)) 
recovers the electron cyclotron resonance fi = (l,i,0) 
for uji = ujc such that (pi = 7r/2. Now suppose that 
62 — k{x + y + z)/3^/^ — ujt, where the square root of 
3 normalizes the unit vector. The explicit solution now 
gives f2 oc (l,e'2'^/3,e-*2'^/3) for 0J2 = i^cji^''^ such that 
02 = arcsin3~^/^ and O2 — t:/A. One can check that 
Eqn. (|55|) reproduces these phase factors using the ap- 
propriate angles for the propagation vector. 

The trigonometric functions extend analytically to the 
complex plane (p ^ (p such that cos^ (p + sin^ (p = 1 is 
preserved. In particular, one can extend the expres- 
sion si'a(p± — ±uj/ujc to the domain w > Wc such that 
cos(p± = i{uj'^ /loI — 1)^/^. The evolution of the complex 
(p+ as a function of the frequency ratio lo/oJc is depicted 
in Fig. [51 as are its sine and cosine. The angle (j)+ pro- 
ceeds along the real axis from to 7r/2 as w goes from 
to We, and then it acquires an imaginary component 
when Ld > LOc- The sine is entirely real, while the cosine 
goes from being real to imaginary as the frequency ra- 
tio exceeds unity; consequently, the unit vector k is itself 
complex in this regime. Nonetheless, it maintains its nor- 
malization k • k = 1 and its orthogonality to the vector 
phasor k • f = 0. We will come back to the interpre- 
tation of k after we solve the conductive d'Alembertian 
equation for its magnitude k. 
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FIG. 2. (Color online.) Shown are the real □ and imaginary 
components of </>+ as a function of the frequency ratio tj/cjc 
in (a) and (b), as well as its sine in (c) and its cosine in (d). 

What remains is to solve Eqn. pi)) for fc(cj), which is 
a standard problem in electrodynamics Having writ- 
ten the propagation vector in terms of its magnitude 
and direction, the d'Alembertian operator simplifies con- 
siderably, — = — u'^/cq. Rearranging factors af- 
ter writing Vg in terms of A yields the scalar relation 
— = ipLQ<juj for DC conductivity a = e^ne/meVei, 
which has the positive root 

k = {uj/co){l + i(j/eo^Y'^ (56) 

with magnitude k — (cj/co)(l + CT^/epW^)^/'* and phase 
5k ~ arctan(tT/eow)/2. When uj < lOc the real propa- 
gation and decay vectors k_|_ and k_ point in the same 
direction such that k — |k+| -|- i|k_|, but when uj > tOc 
the situation gets more complicated because of the com- 

plex form of k = k+ -I- ik_, where |k_|_| = oj^/cj^ and 
|k_| = w^/w^ — 1. In this notation, k+ ^ k+/|k+| as 
each expression refers to a different object; k+ is the real 
part of k, whereas k+/|k+| gives the direction of prop- 
agation according to the real part of k. The direction 
of decay is given by k_/|k_|, which is not aligned with 
the direction of propagation when the frequency ratio 
exceeds unity. Recalling k = k+ + ik_ such that 

exp'(k • r) = exp(-k_ • r) exp*(k+ • r) , (57) 

the complex propagation vector describes the dispersion 
of electromagnetic radiation in a conductive medium with 
skin depth A^- = l/|k_|, wavelength A^^ — 27r/|k+|, phase 
speed Cu: = CLi/|k+|, and group speed — 9a;/9|k+|. 



IV. EVALUATION OF THE NONLINEAR SOLUTION 

The evaluation of the solution to the system of equa- 
tions as a function of the frequency ratio w/wc divides 
naturally into that for the vector phasor described by f 
and for the scalar phasor described by k. The solution 
is translationally invariant, as it is expressed in terms of 
its magnitude V at some arbitrary location identified as 
the origin of the coordinate system. The expression for f 
in Eqn. (|55p depends only upon the frequency ratio and 
the arbitrary azimuth 9^, whereas the expression for k 
depends additionally upon the conductivity a. Let us 
now look at each in turn. 



A. Evaluation of the vector phasor 

The vector phasor f = f+ + if- describes the direc- 
tion of the electron flow at two times separated by one 
quarter of a cycle. Under the conditions of normaliza- 
tion and orthogonality it can be reduced to two scalar 
degrees of freedom corresponding to two of the compo- 
nents of Eqn. p8|) . but it is instructive to evaluate all 
of its components explicitly. The third component of 
Eqn. determines the direction k through its depen- 
dence on (j)± as a function of the frequency ratio uj/uJc- 

Shown in Fig. [3] are the {x, y, z) components of f+ and 
f _ as a function of uj /uJc for various azimuths 6k spanning 
one half revolution and selecting 0+ . The components of 
k can be inferred from Fig. [5J Focusing on the case of 
6k — such that k has no y component, the flow is 
primarily in the y-z plane for w <C Wc, in the x-y plane 
for UJ = UJc, and in the x-z plane for uj ^ uJc- Despite 
having flow in the same plane as propagation when the 
frequency ratio exceeds unity, one can verify that the 
divergence constraint k ■ f = is satisfied. 



B. Evaluation of the scalar phasor 

The scalar phasor exp*((5) for S = k-r— describes the 
behavior of the time harmonic solution in space through 
the complex propagation vector k as a function of the fre- 
quency ratio uj/uJc and propagation azimuth 9k, as well 
as the plasma conductivity a. In the following let us set 
6k = so that propagation and decay are solely in the x-z 
plane, and let us select the (f>^ solution. The electromag- 
netic sector reduces to the scalar equation — = i^qguj 
which determines k according to Eqn. (j56p . thus we need 
to consider what is a reasonable range for the parameter 
a in units of Siemens per meter. 

The conductivity is proportional to the ratio of den- 
sity to collision frequency a oc rie/vei, where v^i de- 
pends upon both density and temperature. The range 
of density and temperature spanned by matter in the 
plasma state covers many orders of magnitude for each. 
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FIG. 3. (Color online.) The vector phasor f, here for <l)+, 
depends upon the frequency ratio uj/lJc and the propagation 
azimuth 8k, equal to in (a) and (b), n/4 in (c) and (d), n/2 
in (e) and (f), and 37r/4 in (g) and (h). On the left is its real 
part f+ , and on the right is its imaginary part f_ . The {x,y, z) 
components are displayed as □, Oi s-nd respectively. 



A common parametrization of the n-T plane, where T 
is the thermal energy, is given by the (electron) Debye 
length \d = (eoTe/riee^)^/^ and Debye density ud = 
(47r/3)A'^ne- Their contours are shown in Fig. 2] over 
a reasonable range of the n-T plane. Also shown are 
contours for the conductivity evaluated as follows. For a 
gaseous plasmal^ one can estimate the collision rate from 
the expression 



1^-' = 6e2(6^m,T3)i/2/(n,e4logA) , 



(58) 



where A = 127r[{eoTe / e"^)^ /rie]^^^ . One can see that the 
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FIG. 4. (Color online.) Shown are contours for conductivity 
(7 (solid), Debye length A_d (dashed), and Debye density no 
(dash-dot) over the n-T plane with temperature expressed in 
units of electron-volts. The contours for a are labeled verti- 
cally, for \o are horizontal, and for no are to the outside of 
the plot. 



conductivity a has only a weak dependence on density 
Ue through log A. A reasonable range for a can thus be 
estimated as 10" < cr < 10^° which spans materials from 
poorly to highly conductive. 

The expression for k depends explicitly on the ratio 
of conduction to displacement current p = a/eQUj = 
tan(2(5fc). In the presence of a background magnetic 
field Bq with magnitude Bq in units of tesla, that ra- 
tio can be rewritten in terms of the cyclotron frequency 
as p = {u!c/uj)po, where 



PO = cr/ eoUJc 



(59) 



is a unitless parameter depending on the conductivity 
and the magnitude of the background field. The value of 
Pq is ultimately what determines the dispersion of elec- 
tromagnetic radiation in a conductive medium subject to 
an external magnetic field. 

To determine the real propagation and decay vectors 
k_|_ and k_ , we construct the complex propagation vector 
k according to its magnitude k and direction k and then 
decompose its real and imaginary parts, k+ -f ik_ = fck. 
The directions for propagation and decay are given by 
the unit vectors k+/|k_|_| and k_/|k_|, which for uj > ojc 
are found not to point in the same direction as a phase in 
addition to Sk appears when k is complex. When uj < uj^ 
the unit vector k is real, and for all cj one finds k • (f+ x 
f_) = l. 

For various po G [10 ^, 10^], the directions of propaga- 
tion and decay as a function of the frequency ratio are 
shown in Fig. \5\ When w < Wc they are both equal to 
k = f+ X f_, which does not depend on po- For uj ^ ujc 
they are nearly but not quite orthogonal, as k^ • k_ > 0, 
and the extent of the transition region does depend on 
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Pq. Attempting to describe the behavior, as oj goes from 
to ujc the propagation and decay vectors both point 
along the normal to the flow plane, which in this case 
9k = goes from x to z in the x-z plane. As lo increases 
beyond the cyclotron frequency, the propagation direc- 
tion first acquires a component along — x before swinging 
back to z. The direction of decay proceeds back towards 
X, and the normal to the flow plane heads towards y in 
the y-z plane. As the ratio po increases, the transition 
region for the propagation direction grows both in terms 
of the swing towards — x and how high a frequency ratio 
is needed before pointing along z, and similarly for the 
direction of decay. 

The magnitudes of the propagation and decay vectors 
determine the wavelength and skin depth as a function of 
the frequency ratio, Ai^ = 27r/|k+| and Ao- — l/|k_|. The 
dispersion relation yields the phase speed c^^ = w/|k+|, 
and its derivative the group speed w,^ = duj/d\\i+\. For 
all uj one can decompose the complex propagation vector 
k as 



kk = (fc+ + ik^)(k+ + ik_) 
= (fc+k+ - fc-k-) + i{k^k^ 



fc+k,) 



(60a) 
(60b) 



and when k_ = one has |k+| = fc+ and |k_| = k^. 
In Fig. [S] we display A^i and Ao- in units of the vacuum 
wavelength Ao — co/(w/27r) as well as c^j and in units 

I 



of the vacuum light speed cq, for po S [10~^, 10^]. 

The transition between conductive and non-conductive 
behavior occurs for a value of po J; 1, such that A^ > A^ 
for a range of frequencies w < Wc. In other words, sig- 
nificant propagation of the oscillation over several wave- 
lengths occurs only for a span of frequencies below the 
cyclotron resonance when the conductivity is sufficiently 
low that dissipation does not destroy the waveform. For 
higher conductivities the skin depth is but a fraction 
of the wavelength, indicating that the amplitude decays 
to almost nothing before a single cycle is realized. For 
UJ < iOc the phase speed is less than the group speed, and 
for UJ > Uc it is greater. Interestingly, for a poor conduc- 
tor Po ^ 1 there exists a range of frequencies below ujc for 
which the group speed exceeds the vacuum light speed. 
Such a case is not unheard of in plasma physics,^- and we 
stress that this result obtains directly from the accepted 
solution to the inhomogeneous d'Alembertian equation 
for conductorsji^ One can verify that the group speeds 
displayed, evaluated from numerical gradients, agree with 
the analytic expression 



2(p2 + 1)3/4^^ 



p sin Sk + (p^ + 2) cos Sk 



(61) 



for frequencies below the cyclotron resonance, and above 
the resonance 



2(p^ + 1)3/4 [(^/u,,)2_sm^SkY^\o 
p cos Sk sin Sk — (p^ -I- 2) sin^ Sk + (aj/wc)^(3p^ + 4) 



(62) 



recalling p = tan(2(5fc). Equation (j6T|) can be derived 
from textbook electrodynamics, while Eqn. is a con- 
tribution of this investigation taking into account k_ . 



V. DISCUSSION 

In discussing the significance of these results, let us be- 
gin by enumerating some of the many effects neglected by 
this analysis. By using the vacuum values for permittiv- 
ity and permeability eo and po i one has ignored both the 
electric and magnetic atomic polarizability of the ions. 
Obviously these effects will have an impact on the eval- 
uation of the wave number k and all quantities derived 
therefrom. For that reason we are not overly concerned 
by the appearance of a group speed in excess of the vac- 
uum light speed at this stage in the development of the 
theory. When the atomic polarizabilities can by absorbed 
into scalar relative permittivity and permeability e and 
p, one simply replaces the vacuum expressions with those 
values; however, in a magnetized plasma, those quanti- 
ties are usually defined as gyrotropic tensors derived from 
the equations of motion. 



I 

The motion of the ions is neglected throughout these 
derivations. To the extent that the ion and electron mo- 
tions decouple, one can repeat the analysis for J = eniVi 
over a range of frequencies around the ion cyclotron fre- 
quency to yield a "slow" partner to the "fast" mode de- 
rived above for the electrons. As the ions (are assumed 
to) carry positive charge, their motion couples to electro- 
magnetic oscillations of opposite circular polarization to 
that for electrons The most interesting case of course is 
when both electrons and ions are allowed to flow. In that 
situation their combined flow can be decomposed into 
net current and momentum oscillations J and K coupled 
through the equations of motion. One expects that low 
frequencies will excite mostly momentum waves, whereas 
high frequencies should drive mostly current waves, with 
separate transition frequencies for positive and negative 
helical polarization. 

An important effect neglected here is the possibility 
of temporal and spatial variation to the electron density, 
as well as the presence of sources or sinks of current, 
(i9t -t- Ve • V)ne = He — (V ■ Ve)ne. It is the effect of a net 
charge density, after all, which pushes current around a 
source driven wire; Hernandes and Assis-^"^ give an exam- 
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FIG. 5. (Color online.) Shown are the unit vectors giving the 
directions of propagation k+ and decay k_ as a function of the 
frequency ratio for various values of po indicated to the lower 
left of each panel. The {x, y, z) components are displayed as 
□ , Oi respectively. 



FIG. 6. (Color online.) Shown on the left are the wavelength 
\u) as □ and the skin depth Ao- as in units of the vacuum 
wavelength Ao as a function of the frequency ratio for various 
Po indicated to the upper right of each panel. On the right 
are the phase speed as □ and the group speed v^j as in 
units of the vacuum light speed Co. 



pie of its derivation in the steady state and Jefimenko^i 
gives an example of its measurement. In Lorenz gauge, 
the net charge density j couples to the scalar potential a 
through the same d'Alembertian operator — D^a = ^qj. 
With respect to the derivation of the (electron) plasma 
frequency lOp = riee^metQ, a non- vanishing j should allow 
for the coupling of the electron fluid to electromagnetic 
oscillations of linear polarization in addition to the circu- 
lar polarization described here. In this context, it might 
be beneficial to decompose the statement of net momen- 
tum conservation Eqn. (I34[) in terms of the Maxwell stress 



tenso r^^'^^ so that the components of the macroscopic 
Lorentz force can be identified. 

As the only external potential apparent in the the- 
ory is that of the background field Bq, what we have 
derived is essentially a model for resonance within the 
electron fluid, not its response to a driving potential at 
arbitrary frequency. In order to describe the plasma re- 
sponse to an external device, such as a radio frequency 
antenna, the potential from that device must be calcu- 
lated and included explicitly in the mathematics. For the 
common configuration of a plasma confined by a cylindri- 
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cal vessel ^"-SS, one also must be explicit with the transi- 
tion from conductive to dielectric material at the bound- 
ary such that rie — 0. Nonetheless, the model derived 
here is based upon the same simplifying approximations 
as commonly used in the theoretical description of such 
devices^'^ii^'^ and sheds some light on the electromag- 
netic behavior of ionized material. 

What we have learned is that one need not neglect the 
nonlinear term in the equation of motion to find an an- 
alytic solution when the ion contribution is included as 
a constraint relating the electron flow to the potential it 
generates. Working with the electromagnetic potential, 
neither is Gauss's law neglected; rather, it is the possible 
occurrence of a net charge density which is neglected in 
the neutral fluid approximation. Our main result is that 
propagation occurs in an ionized medium under the given 
conditions only for a range of frequencies below the elec- 
tron cyclotron resonance when the conductivity is suffi- 
ciently low that the displacement current outweighs the 
conduction current. For a highly conductive medium, the 
skin depth dominates the wavelength such that the am- 
plitude of oscillation is attenuated. Below the cyclotron 
resonance the propagation and decay vectors point along 
the normal to the flow plane, whereas above the reso- 
nance they point in different directions while still respect- 
ing continuity. 



VI. CONCLUSION 

In closing, we hope that this article has contributed 
to the understanding of the analysis of dispersion in 
an ionized medium subject to a background magnetic 
field. The solution presented here satisfies the classical 
Maxwell equations, which may be expressed succinctly in 
terms of geometric formsi^ii^ as d ★ dA = J for ddA = 0, 
as well as both the ion and electron equations of motion 
under the given simplifying approximations. From the 
plane wave expansion one should be able to construct so- 
lutions of arbitrary geometry. The nonlinear coupling of 
the electron flow to its own potential leads to a model 
which specifles a relation between the frequency of os- 
cillation and the inclination of the propagation vector 
from the plane normal to the background held. The cir- 
culatory flow given by the vector phasor describes the 
magnetization of the electron fluid in the presence of a 
propagating electromagnetic oscillation. 

To describe systems of physical interest, this model 
must be extended to a multispecies formalism which 
allows the ions to flow, incorporates the effects of in- 
traspecies collisions through the thermal stress tensors, 
and accounts for a net charge density. Nonetheless, this 
simplified model reveals some interesting physics, namely 
that the propagation of electromagnetic radiation in a 
uniform plasma subject to a constant magnetic field in- 
duces a rotation of the magnetization away from the axis 
of the background field. The ratio of oscillation and 
cyclotron frequencies determines the inclination of the 



propagation vector, which is orthogonal to the plane of 
flow when that ratio is less than unity. The solution 
decays in space according to the ratio of free to displace- 
ment current, which depends upon the conductivity of 
the medium and the frequency of oscillation. 

The motivation for this investigation is the question of 
whether helicon waves really exist, to which our answer 
is yes, they do exist, just not the ones given by previous 
derivations. Working in the potential formulation simpli- 
fles the analysis by reducing the electromagnetic sector to 
its physical degrees of freedom, so that the nonlinear self- 
interaction between the plasma current and the potential 
it generates can be addressed in the equation of motion. 
Under the stated conditions, an analytic solution can be 
found to the nonlinear system of equations describing res- 
onance within the electron fluid of a uniform, stationary 
plasma subject to a constant background magnetic fleld. 
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